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Abstract 

We present preliminary results of a new measurement of the b quark fragmentation 
function in Z° decays using a novel kinematic B hadron energy reconstruction technique. 
The measurement is performed using 150,000 hadronic Z° events recorded in the SLD 
experiment at SLAC between 1996 and 1997. The small and stable SLC beam spot and 
the CCD-based vertex detector are used to reconstruct topological I?-decay vertices with 
high efficiency and purity, and to provide precise measurements of the kinematic quantities 
used in this technique. We measure the B energy with good efficiency and resolution over 
the full kinematic range. We compare the scaled B hadron energy distribution with several 
functional forms of the B hadron energy distribution and predictions of several models of 
h quark fragmentation. Several functions including JETSET + Peterson are excluded by 
the data. The average scaled energy of the weakly decaying B hadron is measured to be 
x B = 0.719 ± 0.005 (stat) ± 0.007 (syst) ± 0.001 (model) (preliminary). 
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1. Introduction 



The fragmentation process which transforms colored partons into colorless hadrons is 
typically characterized by the fragmentation function. The b quark fragmentation is 
of special importance in the study of quark fragmentation because the large b quark 
mass provides a natural mass scale in QCD calculations and allows the application of 
perturbative QCD and heavy quark expansion, and can help to extract non-perturbative 
effects in fragmentation, which is the least understood part of the fragmentation process. 

According to the factorization theorem, the heavy quark fragmentation function can 
be described as a convolution of perturbative and non-perturbative effects. For the b 
quark, the perturbative calculation is in principle understood |I|, 0, [| f|, Nonpertur- 
bative effects have been parametrized in both model-dependent || |7|, [8|, ||, [10], [l| and 
model-independent approaches J!], |], The fact that several models are yet to be 
experimentally tested is an indication of a lack of precise and conclusive experimental 
results, if not theoretical understanding. 

It is indeed experimentally challenging to measure the b quark fragmentation function 
to a level of precision sufficient to distinguish among the various models. Since the b quark 
fragmentation function is the probability distribution of the fraction of the momentum of 
the b quark carried by the B hadron, the most sensitive experimental determination of 
the shape of the b fragmentation function is expected to come from a precise direct mea- 
surement of the B hadron energy (or momentum) distribution. The difficulty in precisely 
measuring the B hadron energy distribution stems mostly from the fact that most of the 
B decays can only be partially reconstructed, causing a significant fraction of the B energy 
to be missing from the B decay vertex. Recent direct measurements at LEP [11, 12] and 



SLD JT3| have used overall energy-momentum constraints and calorimetric information to 
extract this missing energy in a sample of semi-leptonic B decays. These measurements 
suffer from low statistics as well as poor B energy resolution at low energy, and hence 
have a relatively weak discriminating power between different shapes of the fragmentation 
function. Indirect measurements such as the measurements of the lepton spectrum 
and charged multiplicity have been used to constrain the average B energy. However, 
these measurements are not very sensitive to the shape of the energy distribution. 

Here we report preliminary results of SLD's new measurement of the B hadron energy 
distribution. We developed a novel kinematic technique which uses only charged tracks 
associated with the B vertex and the B flight direction to reconstruct individual B hadron 
energy with good resolution over the full kinematic range while achieving an efficiency 
much higher than previous measurements. 



2. B hadron Selection 



A general description of the SLD detector can be found elsewhere [fT6| , |17|1 The excellent 
tracking and vertexing capabilities at SLD [[TSJ are exploited in the reconstruction of B 
decays in Z° —>■ b b events. 

A set of cuts is applied to select hadronic Z° events well-contained within the detector 
acceptance. The efficiency for selecting a well-contained Z° — > qq(g) event is estimated to 
be above 96% independent of quark flavor. The selected sample comprise 111,569 events, 
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with an estimated 0.10±0.05% background contribution dominated by Z° — > r + r~ events. 

The B sample is selected using a topological vertexing technique based on the detection 
and measurement of charged tracks, which is described in full detail in Ref. [HI. The 



topological vertexing algorithm [|19|] is applied separately to the set of "quality" tracks in 
each hemisphere (defined with respect to the event thrust axis). 

When a candidate vertex is found, tracks not associated with this "seed" vertex are 
attached to the vertex if they are more likely to have originated from this vertex than from 
the IP. This track-attachment procedure is tuned to minimize false track-vertex associ- 
ations to the vertex. Attaching a false track to the vertex affects the vertex-kinematics 
more than failing to associate a genuine track originated from the vertex, and hence can 
cause significant degradation in the reconstructed B energy resolution. On average, this 
procedure attaches 0.8 tracks to each seed vertex; about 92% of the reconstructed tracks 
which originated from the 5-decay are associated with the reconstructed vertex, and 98% 
of the vertex-associated tracks are true B decays tracks. 

The mass of the reconstructed vertex, 
M c h, is calculated by assigning each track the 
charged-pion mass. Because of the tiny SLC 
IP error and the excellent vertex resolution, 
the B flight direction pointing along the line 
joining the IP and the secondary vertex is 
well-measured. Therefore the transverse mo- 
mentum P t of tracks associated with the ver- 
tex relative to the B flight direction is also 
well-measured. The mass of the missing par- 
ticles can then be partially compensated by 
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Figure 1 Distribution of the recon- 
structed Pi-corrected vertex mass in 
the 1996-97 data (points). Also shown 
is the prediction of the Monte Carlo 
simulation, for which the flavor com- 
position is indicated. 



using P t 

M Pt = 



to 



M 2 ch 



form the u P t corrected mass", 
To minimize ef- 



P? + 



Pt 



fects of large fluctuations of P t at short de- 
cay length, the minimum transverse momen- 
tum (which is varied within the la limits con- 
straining the axis at the measured interaction 
point (IP) and reconstructed seed vertex) is 
calculated in order to determine Mp t . Figure 
1 shows the distribution of the P t -corrected 
mass (points) for the 32,492 accepted hemi- 
spheres in the data sample, and the corresponding simulated distribution. To obtain a 
high purity B sample, B hadron candidates are selected by requiring Mp t > 2.0 GeV/c 2 . 
A total of 19,604 hemispheres are selected, with an estimated efficiency for selecting a true 
P-hemisphere of 40.1%, and a sample purity of 98.2%. The contributions from light-flavor 
events in the sample are 0.15% for primary u,d and s events and 1.6% for c events. 



3. B Energy Reconstruction 

Since the sum of the charged track energy at the B vertex, E c h, is known, we are only 
concerned with finding the energy of particles missing from the B vertex. 
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Figure 2. The relative deviation of the max- 
imum missing mass from the true missing 
mass for Monte Carlo simulated B hadron 
decays, which is divided into three cat- 
egories: B° and B^ (open), B® (cross- 
hatched), and Afe (dark hatched). 
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Figure 3. Distribution of the reconstructed 
M^ max for the selected vertices in the 1996- 
97 data (points). Also shown is the predic- 
tion of the Monte Carlo simulation. 



Given reconstructed B vertex, an upper bound on the mass of the missing particles 



from the vertex is found to be MQ max 



Ml 



2M B ^M 2 ch + P t 2 + M c \, where we assume 
the true mass of the B hadron decayed at the vertex, Mb, equals the B° meson mass. 
Since the true missing mass M^ ue is often rather close to M 0max (Figure 2), M 0max is 
subsequently used as an estimate of Mq™ 6 (Mq U£ is set to if the reconstructed M 0rnax is 
negative) to solve for the longitudinal momentum of the missing particles from kinematics: 
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Mb_ 



--{Ml 



2(MX+Pf) 



PZlp 

1 I 



chl, 



and hence the missing B energy from the vertex, E . The B hadron energy is then 
E B = E ch + E . Since < M l ™ e < M 0max , the B energy is well-constrained when M 0max 
is small. In addition, most uds and c backgrounds are concentrated at large M 0max . We 
choose an ad hoc upper cut on the M^ max to achieve a nearly x^-independent B selection 
efficiency. Figure 3 shows the distribution of Ml max after these cuts, where the data and 
Monte Carlo simulation are in good agreement. A total of 1920 vertices in the 1996-97 
data satisfy all selection cuts. Figure 4 shows the distribution of the reconstructed scaled 
weakly decaying B hadron energy for data and Monte Carlo. The overall B selection 
efficiency is 3.9% and the estimated purity is about 99.5%. The efficiency as a function 
of x l g ue is shown in Figure 5. We examine the normalized difference between the true and 
reconstructed B hadron energies for Monte Carlo events. The distribution is fitted by a 
double Gaussian, resulting in a core width (the width of the narrower Gaussian) of 10.4% 
and a tail width (the width of the wider Gaussian) of 23.6% with a core fraction of 83%. 
Figure 6 shows the core and tail widths as a function of x% ue . The core width depends 
only weakly on the true x B , another feature that makes this method unique. 
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Figure 4. Distribution of the reconstructed scaled B hadron energy for data(points) and 
simulation (histogram). The solid histogram shows the non-66 background. 
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Figure 5. Distribution of the efficiency as a Figure 6. The fitted core and tail widths of 



function of the true B energy. 



the B energy resolution function of the 
true scaled B hadron energy. 
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4. Tests of Functional Forms and Models 



After background subtraction, the distribution of the reconstructed scaled B hadron en- 
ergy is compared with a set of ad hoc functional forms of the xb distribution in order 
to estimate the variation in the shape of the xb distribution. For each functional form, 
the default SLD Monte Carlo is re-weighted and then compared with the data bin-by-bin 
and a \ 2 is computed. The minimum \ 2 is found by varying the input parameter (s). The 
Peterson function ||, two ad hoc generalizations of the Peterson function |i2|| (ALEPH 1 
and 2) and a 7th-order polynomial [] are consistent with the data. We exclude the func- 
tional forms described in BCFY Q, Collins and Spiller [10|, Kartvelishvili ||, Lund §J 
and a power function of the form f(x) = x a (l — x 
in Table 1 and 2. 



,/3 



The result is shown in Figure 6 and 
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Figure 7. Each figure shows the background-subtracted distribution of reconstructed B 
hadron energy for the data (points) and for the simulation (histograms) based on the 
respective optimised input fragmentation function. The x 2 fit uses data in the bins between 
the two arrows. 



'The behavior of this polynomial is rather unphysical at low Xb and will not be considered hereafter 
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Function D(x) Reference 

ALEPH 1 ^^(l-i-T^r 2 
ALEPH 2 1(1 - § - 



BCFY [1 - (1 Jffaf [3 + ^=x(-xYfM I 

Collins and Spiller (i^ + + x 2 )(l - 1 - y^)" 2 

Kartvelishvili ei a/. x Qfc (l — x) 

Lund ^(l — x) a exp(— bm\/x) @ 

Peterson et oZ. |(1 - | - jt^) -2 [| 

Polynomial x(l — x){l + Sf =1 PjX l ) (see text) 

Power x a (l — x)P (see text) 



Table 1. Fragmentation functional forms used in comparison with the data. For the BCFY 
function /i(r) = 3(3 - 4r), f 2 (r) = 12 - 23r + 26r 2 , f 3 (r) = (1 - r)(9 - llr + 12r 2 ), 
and /i(r) = 3(1 — r) 2 (l — r + r 2 ). A polynomial function and a power function are also 
included. 



Function 



XVdof 



Parameters 



ALEPH 1 
ALEPH 2 

BCFY 

Collins and Spiller 
Kartvelishvili et al. 
Lund 

Peterson et al. 
Polynomial 



Power 



15.2/15 
23.7/15 

52.3/16 
54.3/16 
79.6/16 
139.1/15 

26.0/16 
14.4/12 



78.5/15 



0.860 



+o.uiy 

-0.018 



(x B ) 



d = 0.019 ±0.002 



0.938 



+0.039 
-0.034 



d = 0.036 ±0.002 



0.2316 



+0.0092 
-0.0088 



£6 



n n44+ - 005 

U.U++_ 004 



a b = 4.15 ±0.11 



2.116 



+0.118 
-0.114 



bm\ = 0.408ig:8?g 
Q> = 0.0338t8:oo22 
p x = -12.4 ±0.4 

p 2 = 58.7 ± 1.9 
p 3 = -130.5 ±4.2 
p 4 = 136.8 ±4.3 
p 5 = -53.7 ± 1.8 
a = 3.91181 

P = o.894i8: 8? 



0.718±0.005 
0.720±0.005 

0.713±0.005 
0.714±0.005 
0.720±0.004 
0.720±0.005 

0.719±0.005 

(see text) 



0.722±0.005 



Table 2. Results of the \ 2 fit of fragmentation functions to the reconstructed B hadron 
energy distribution after background subtraction. Minimum x 2 , number of degrees of 
freedom and coresponding parameter values are listed. Errors are statistical only. 
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We then test several heavy quark fragmentation models. Since the fragmentation 
functions are usually functions of an experimentally inaccessible variable (e.g. z = 



(E+p\\)h/ (E+p\\)q), it is necessary to use a Monte Carlo generator such as JETSET [^C 
to generate events according to a given input heavy quark fragmentation function. The 
resulting B energy distribution is then used to re-weigh the Monte Carlo distribution 
before comparing with the data. The minimum \ 2 is found by varying the input param- 
eter (s). Within the context of the JETSET Monte Carlo, Bowler and the Lund §| 
models are consistent with the data, while Peterson || model is found to be inconsistent 
with the data. 



5. Systematic Errors 

We have considered both detector and physics modelling systematics. The dominant 
systematic error is related to charged track transverse momentum resolution smearing, 
which has been evaluated conservatively and can be reduced with a detailed study. All 
physics systematics are rather small. Other relevant systematic effects such as by varying 
the event selection cuts and the assumed B hadron mass are also found to be small. 
In each case, conclusions about the shape of the B energy distribution hold, and the 
systematics in the average B hadron energy is added in quadrature to obtain the total 
systematics. 



6. Conclusion 

Taking advantage of SLC's small beam-spot and SLD's high vertex resolution, we have 
developed a new kinematic technique to measure, for the first time, the B hadron energy 
distribution in Z° decays with good resolution over the full kinematic range. Using 1996-97 
data, we exclude several functional forms of the B energy distribution and the JETSET + 
Peterson fragmentation model. The mean of the scaled weakly decaying B hadron energy 
distribution is obtained by taking the average of the means of the three functional forms 
which are found to be consistent with the data. The r.m.s. of the three means is regarded 
as a minimum error on model-dependence. We find 

< x B >= 0.719 ± 0.005{stat) ± 0.007{syst) ± O.OOl(modeZ) 

where the small model-dependence error indicates that < xb > is relatively insensitive 
to the allowed forms of the shape of the fragmentation function. The precision in the 
measured average B hadron energy represents a substantial improvement over previous 
direct measurements. All results are preliminary. 
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